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A I  computational  study  has  been  performed  to  determine  the  aerodynamics 
of  a  cluster  of  roimd  parachutes  using  computational  fluid  dynamics 
(GFD).  The  results  given  here  are  the  first  predictions  of  descent 
cliaracteristics  for  a  cluster  of  three  half-scale  C-9  parachutes.  In 
p^icular,  the  results  include  the  aerodynamic  flow  field  and  geometry  of 
thb  parachute  cluster  with  assumed  shapes  for  the  individual  canopies. 
Computed  results  have  been  obtained  using  both  structured  and 
unstructured  numerical  techniques.  The  computed  pressure  over  the  inner 
and  outer  surfaces  of  a  single  canopy  in  the  cluster  is  used  to  calculate  the 
net  forces  and  moments  acting  on  the  canopy.  A  manual  iterative  process 
isj  used  to  determine  the  expected  stable  configuration  for  the  cluster 
geometry  by  determining  the  condition  at  which  the  net  moment  about  the 
payload  (origin)  is  zero.  The  corresponding  angular  location  of  the 
canopies  predicted  by  the  CFD  computations  is  compared  with  the 
available  experimental  results  and  is  found  to  agree  well  with  the  data.  It  is 
shown  that  significant  progress  has  been  made  in  determining  the  terminal 
descent  flow  field  characteristics  of  a  particular  parachute  cluster 
configuration.  These  computational  solutions  provide,  for  the  first  time, 
ah  understanding  of  the  flow  field  in  and  around  parachute  clusters. 
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COMPUTATIONAL  FLUID  DYNAMICS  MODELING  OF  PARACHUTE  CLUSTERS 


1.  INTRODUCTION 

Parachute  clusters  have  been  used  in  a  wide  range  of  applications  throughout  their  history, 
and  the  process  utilized  for  their  development  and  design  has  significantly  improved  over  the 
years.  Parachute  cluster  systems  were  developed  to  avoid  the  use  of  excessively  large  single 
canopies.  A  practical  limit  for  the  utility  of  single  round  canopies  is  approached,  as  the  users’ 
desire  for  recovering  larger  payload  weights  increases.  In  many  applications,  a  cluster  of 
parachutes  has  a  number  of  benefits  over  a  large  single  canopy.  These  include  the  ability  to  rig 
and  manufacture  the  system.  Clusters  also  have  the  advantage  of  backup  protection.  This  is 
evident  when  a  single  cluster  chute  failure  is  considered.  The  system  imder  a  smaller  subset  of 
the  canopies  still  has  a  lower  impact  velocity  than  a  free  fall.  Clustered  parachute  systems  have 
excellent  stability  characteristics  compared  to  most  single-canopy  systems.  This  is  partly 
because  of  the  terminal  descent  shape  of  multiple  bluff  bodies  versus  a  single  bluff  body.  Of 
course,  there  are  some  undesirable  features  of  parachute  clusters,  which  include  the  difficulty  in 
obtaining  a  time-sequenced  opening  of  all  canopies  together.  This  phenomenon  is  known  as  “lead 
lag,”  in  which  one  canopy  opens  quicker  than  the  next  canopy.  This  generally  forces  the 
individual  parachutes  in  the  cluster  to  be  stronger  than  would  be  needed  for  a  consistent  opening 
of  all  canopies  together.  Another  disadvantage  of  parachute  clusters  is  that  the  total  drag  from 
each  canopy  in  the  cluster  is  less  than  the  drag  that  would  be  experienced  by  the  equivalent  single 
canopy  at  the  same  descent  velocity.  This  reduction  in  drag  is  attributable  to  the  shape  about 
which  clustered  parachute  systems  oscillate.  Each  canopy  in  the  cluster  is  most  often  positioned 
at  some  angle  of  attack  to  the  trajectory  path  of  the  system,  which  reduces  the  drag  produced  by 
each  canopy.  The  optimization  of  a  cluster  configuration  can  result  in  a  trade-off  between  the 
overall  height  of  the  system,  which  affects  the  opening  time,  etc.,  to  the  lengths  used  for  the 
risers. 

Most  cluster  applications  described  in  the  literature  are  National  Aeronautics  and  Space 
Administration  (NASA)  systems,  including  the  Apollo  recovery  systems  and  the  space  shuttle 
solid  rocket  booster  recovery  systems.  However,  parachute  clusters  are  employed  for  many 
other  applications,  including  numerous  military  applications.  Military  systems  that  use  clusters 
include  extraction  systems  and  low-velocity  airdrop  systems  (LVADS)^  for  cargo,  which 
consistently  deliver  payloads  as  heavy  as  60,000  pounds  (see  Figure  1).  The  G-1 1  and  G-12  are 
the  most  common  military  parachutes  that  are  used  for  cargo  platforms  in  a  variety  of 
configurations  based  on  the  payload  weight.  For  example,  a  60,000-pound  system  uses  a  cluster 
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of  three  28-foot-diaineteij  extraction  canopies  that  extract  twelve  100-foot-diameter  G-1 1 
canopies.  Lower  weights  use  fewer  canopies  such  as  a  20,000-pound  system,  which  typically 
uses  a  single  28-foot  ring  slot  extraction  canopy  to  extract  four  100-foot-diameter  G-11  canopies. 
One  design  challenge  of  these  systems  is  the  desire  for  commonality.  The  users  want  systems 
that  can  be  employed  for  a  wide  range  of  weights;  they  also  want  to  minimize  the  number  of 
separate  items  needed  in  inventory.  Therefore,  it  is  desirable  to  have  a  predictive  tool  that  can  be 
used  to  accurately  predict  the  terminal  descent  velocity,  configuration,  and  stability  of  a  given 
system  for  a  wide  range  of  payload  weights.  A  computational  tool  that  models  the  terminal 
descent  characteristics  of  a  cluster  of  parachutes  is  a  technology  that  is  needed  by  parachute 
designers  and  engineers.  A  joint  effort  between  the  U.S.  Army  Soldier  Systems  Command, 
Natick  Research  Development  and  Engineering  Center  (Natick),  and  the  U.S.  Army  Research 
Laboratory  (ARL)  to  develop  this  computational  tool  is  now  under  way.  Current  cluster 
parachute  systems  are  overly  designed  and  often  poorly  optimized  because  the  interacting  flow 
field  associated  with  parachute  clusters  is  not  understood. 


As  with  most  parachute  systems,  the  tools  available  to  predict  performance  and  to  design  a 
clustered  parachute  system  are  based  on  empirical  data  and  comparison  with  proven  systems.  A 
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lot  of  experience  helps.  Some  fiandamental  experimental  work  has  been  done  with  clustered 
parachute  systems.^’^  However,  very  little  is  understood  about  the  flow  field  characteristics  in 
and  around  clustered  parachutes.  One  goal  of  this  work  is  to  provide  a  fundamental 
understanding  of  parachute  cluster  flow  fields.  The  ability  to  model  the  opening  process  of  a 
cluster  of  parachutes  is  extremely  complex  and  is  not  a  current  goal  of  this  work.  The  goal  of  this 
work  is  to  model  terminal  descent  characteristics  and  the  stability  of  the  system  once  all  canopies 
are  open.  The  ultimate  goal  is  to  develop  a  predictive  tool  that  can  be  used  to  optimize  a 
parachute  cluster’s  terminal  descent  characteristics.  This  will  include  the  ability  to  predict  the 
effect  of  modifying  riser  lengths  and  suspension  line  lengths  on  the  canopies. 

The  ability  to  accurately  simulate  the  flow  field  around  a  parachute  cluster  is  a  complex 
problem.  The  aerodynamic  characteristics  associated  with  a  single  or  cluster  of  parachutes  in  the 
terminal  descent  phase  is  extremely  complex  to  model.  The  complexity  arises  largely  from  the 
fact  that  the  flow  field  depends  on  the  canopy  shape,  which  itself  depends  on  the  flow  field.  A 
correct  model  must  include  the  coupled  behavior  of  the  parachute  system’s  structural  dynamics 
with  the  aerodynamics  of  the  surrounding  flow  field.  A  coupled  model  is  required  to  determine 
the  terminal  descent  characteristics  of  clustered  parachutes  including  velocity,  shape,  drag, 
pressure  distribution,  and  the  other  flow  field  characteristics.  As  a  starting  point,  the  present 
research  has  focused  on  the  use  of  computational  fluid  dynamics  (CFD)  to  gain  a  basic 
understanding  of  the  interference  aerodynamic  flow  fields  associated  with  parachute  clusters. 
Prior  work  included  CFD  modeling  of  single  “axisymmetric”  and  three-dimensional  (3-D) 
canopies.  In  this  case,  CFD  techniques  were  used  to  model  a  fixed  shape  single  canopy  in  3-D. 
The  results  were  manually  coupled  to  a  static  structural  code  that  predicted  the  canopy  shape 
based  on  a  CFD-supplied  canopy  surface  pressure  distribution.^®’*^  These  solutions  also 
provided  a  separately  produced  set  of  flow  field  predictions  from  a  completely  different  CFD 
software  tool.  Computed  results  were  compared  with  numerical  results,  which  were  obtained 
earlier  and  used  with  other  concurrently  developed  single-canopy  opening  models.  For  CFD 
modeling  of  parachute  clusters,  we  have  chosen  to  model  a  cluster  of  three  half  scale  C-9  flat 
circular  parachutes.  The  rationale  for  choosing  this  configuration  was  our  ability  to  compare 
some  aspects  of  the  solutions  to  recently  completed  experimental  studies  of  this  cluster 
configuration  during  controlled  conditions  by  Lee,  Lanza,  and  Buckley.*  Once  the  model  is 
validated  against  these  experiments,  the  procedure  can  be  applied  to  other  cluster  systems.  This 
report  describes  the  approach  chosen  to  model  the  cluster  configuration,  the  gridding  procedure 
required,  the  solution  scheme,  the  method  of  comparison  to  the  experimental  work,  and  planned 
future  direction. 
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2.  SOLUTION  TECHNIQUE 


As  with  our  single*canopy  models,  it  is  assumed  that  the  shape  and  configuration  of  the 
cluster  canopies  are  known.  The  initial  surface  shapes  used  for  each  canopy  (axisymmetric  about 
each  canopy  in  the  currdnt  model)  are  the  steady  state  shapes  determined  fi'om  previous 
solutions  of  the  single-cmopy  model  falling  at  approximately  (within  ±1  fi/sec)  the  same  terminal 
descent  velocity,  The  symmetric  shape  of  each  canopy  is  an  approximation,  and  the  accuracy 
of  this  assumption  will  he  studied  by  using  a  3-D  structural  finite  element  program*^  by  feeding  a 
fiill  3-D  representation  of  the  canopy  with  the  CFD  predicted  surface  pressure  distribution  in  the 
near  future.  It  is  anticipated  that  a  few  iterations  will  be  necessary  to  converge  on  the  correct 
canopy  shape  as  was  done  previously  with  the  single-canopy  models.^  ^ 


A  canopy  surface  grid  (for  a  given  cluster  configuration)  at  a  prescribed  polar  angle  (defined 
as  ^)  is  generated  by  u^ing  a  Formula  Translator  (FORTRAN)  code  to  rotate  each  node  point  of 
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the  single  parachute  geometry  in  3-D  about  a  user-defined  confluence  point.  The  polar  angle,  (j) , 
is  measured  from  the  vertical  y-axis  with  the  xz  plane  being  parallel  to  the  ground.  This  fully 
defines  the  shape  of  the  system  being  modeled.  To  modify  suspension  line  lengths,  a  new  single¬ 
canopy  geometry  must  be  generated.  Then  the  riser  lengths  can  be  modified  by  moving  the 
payload  confluence  point  and  rotating  the  geometry  about  the  vertical  axis.  For  the  present 
problem,  four  separate  canopy  surface  orientations  were  considered  and  a  full  3-D  CFD  grid  was 
generated  for  each  case.  For  each  configuration,  the  CFD  code  was  run  with  the  same  boundary 
conditions  and  the  same  user-supplied  in-flow  velocity  (which  is  the  experimentally  determined 
system’s  terminal  descent  velocity).  This  in-flow  velocity  is  defined  to  be  the  velocity  that  was 
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measured  from  the  previously  referenced  experiments.  Solutions  for  each  of  the  four 
configurations  were  obtained,  and  the  canopy  surface  pressure  distributions  were  extracted  from 
the  solution  output.  Note  tliat  the  final  solution  is  a  function  of  the  in-flow  velocity,  which 
ultimately  defines  the  payload  weight  for  the  model. 

The  results  from  each  of  the  configurations  were  used  to  determine  the  correct  configuration 
(polar  angle  for  each  caaopy)  for  a  given  cluster  system.  The  canopy  surface  pressure 
distribution  predicted  by  the  CFD  solution  is  extracted  from  the  overall  3-D  volume  solution.  A 
summation  procedure  is  used  to  calculate  the  net  forces  on  the  fixed  canopy.  The  same  numerical 
fluxes  used  in  the  code  for  the  parachute  surfaces  (inner  and  outer)  are  used  in  the  force  and 
moment  computations.  The  area  of  each  triangle  is  multiplied  by  the  pressure  at  the  cell  center  to 
determine  the  net  force  acting  on  each  triangular  cell.  The  forces  and  moments  are  summed  over 
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all  interior  and  exterior  canopy  triangular  cells.  The  result  of  this  calculation  yields  the  net  forces 
and  moments  acting  on  the  canopy.  The  moment  calculation  is  performed  to  determine  the 
tendency  of  the  cluster  configuration  to  rotate  toward  or  away  from  the  system’s  vertical  y-axis. 
If  the  configuration  consists  of  canopies  that  are  too  close  to  the  vertical  axis,  the  summation  of 
moments  produces  a  net  desire  for  the  configuration  to  separate  from  the  vertical  axis.  A  plot  of 
the  net  moment  for  each  configuration  as  a  function  of  <t>  (the  polar  angle  for  each  canopy  in  the 
cluster)  then  yields  the  predicted  steady  state  configuration  for  the  system  at  the  prescribed  in¬ 
flow  velocity.  This  is  approximated  as  the  value  of  0  at  which  the  net  moment  is  zero.  This 
value  of  0  is  the  computed  steady  state  terminal  descent  polar  angle  about  which  the  system 
oscillates.  Concurrently,  the  net  vertical  drag  produced  by  the  system  is  determined  as  the  net 
vertical  force  of  one  canopy  multiplied  by  the  number  of  canopies.  This  net  force  is  the 
computed  payload  weight  that  will  descend  at  the  terminal  velocity  used  in  the  CFD  calculation. 
Therefore,  by  repeating  this  process  for  a  variety  of  in-flow  velocities,  a  complete  understanding 
of  the  system  for  a  range  of  payload  weights  can  be  determined.  Note  that  the  CFD  solutions  do 
not  reach  a  completely  constant  velocity  and  pressure  field  because  of  vortex  shedding,  etc. 
Therefore,  it  may  be  necessary  to  repeat  these  calculations  over  a  period  of  time  once  the  CFD 
solution  has  evolved. 

2.2  Computational  Fluid  Dynamics  Models 

The  computational  method  used  in  the  initial  analysis  solves  the  incompressible  Navier- 
Stokes  equations  in  3-D^^  generalized  coordinates  for  low-speed  flows.  This  technique  is  based 
on  the  method  of  artificial  compressibility  and  an  upwind  differencing  scheme.  The  pseudo¬ 
compressibility  algorithm  couples  the  pressure  and  velocity  fields  at  the  same  time  level  and 
produces  a  hyperbolic  system  of  equations.  The  upwind  differencing  leads  to  a  more  diagonal 
system  and  does  not  require  a  user-specified  artificial  dissipation.  The  viscous  flux  derivatives 
are  computed  using  central  differencing.  This  code  is  capable  of  computing  both  steady  state  and 
time  accurate  flow  fields.  The  method  is  a  finite  volume  implementation  of  high  accuracy  upwind 
schemes  embedded  in  a  multizone,  structured  grid  framework. 

The  governing  equations  are  numerically  represented  and  solved  using  a  nonfactored  Gauss- 
Seidal  line-relaxation  scheme.  This  maintains  the  stability  and  allows  a  large  pseudo-time  step  to 
be  taken  to  obtain  steady  state  results.  Details  of  the  numerical  method  are  given  in  Pan  and 
Chakravarthy.^^  For  computation  of  turbulent  flows,  a  turbulence  model  must  be  specified.  The 
structured  flow  solver  uses  a  two-equation  k-e  turbulence  model.  In  this  model,  two  trzinsport 
equations  are  solved  for  the  turbulent  kinetic  energy  and  the  turbulent  dissipation  rate.  It  is 
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solved  using  the  same  Gauss-Seidal-type  line-relaxation  scheme  used  to  solve  the  mean  flow 
equations.  ■ 

i 

Another  evolving  technology  that  has  high  potential  in  reducing  the  preprocessing  (grid 
generation)  time  is  unstructured  methodology.  A  recently  developed  unstructured  flow  solver 
by  Rockwell  Science  Center  is  used  in  the  present  parachute  cluster  problem.  It  uses  a 
generalized  Lax-Wendro|lf  scheme  for  solving  Euler  equations  on  unstructured  grids.  This  is  a 

I 

compact  second  order  scheme  that  employs  upwinding  even  in  the  computation  of  the  gradients 
of  the  dependent  variables.  The  generalized  Lax-Wendroff  scheme  is  less  involved,  much  more 
compact,  and  has  the  potential  to  be  a  good  scheme  for  a  wide  range  of  flows.  Also,  the 
compactness  of  the  scheme  enhances  its  stability  and  enables  use  of  larger  values  for  the  time 
step. 

This  approach  employs  a  multilevel  time-stepping  scheme.  The  second  order  scheme  uses 
a  two-time-level  discretization.  For  the  first  fractional  step,  and  are  set  equal  to  the 

corresponding  cell-centejr  values.  Here,  q  contains  all  dependent  flow  variables.  Roe’s 
approximate  Riemann  solver  is  used  to  obtain  a  single  value  for  q(=q  J  at  the  center  of  a  cell. 

For  the  second  fractional  step,  q^  and  qj^  (left  and  right  states)  are  obtained  from  a  Taylor  series 

expansion  of  q  about  the;  centroid  of  the  corresponding  cell.  The  algorithm  just  described  can  be 
considered  as  a  generali^tion  of  Lax-Wendroff  upwind  integration,  since  it  reduces  to  the  Lax- 
Wendroff  scheme  for  uniform  rectangular  hexahedral  cells.  The  solver  does  not  require  any 
information  regarding  the  shapes  of  the  cells  and  the  faces.  This  property  of  the  solver  renders  it 
suitable  for  multigrid  techniques.  In  order  to  be  able  to  compete  with  integral  methods  that  are 
generally  employed  in  the  design  process,  a  multigrid  version  of  the  scheme  described  was 
developed  to  improve  the  convergence  rate.  This  multigrid  version  of  the  code  has  been  used  in 
the  present  parachute  CFD  modeling. 

j 

3,  MODEL  GEOMETRY  AND  COMPUTATIONAL  GRID 

The  CFD  model  was  used  to  determine  the  3-D  terminal  descent  characteristics  of  a  cluster 
of  three  C-9  canopies,  'fhe  constructed  diameter  is  14  feet  and  suspension  lines  are  12  feet  long. 
The  single  half-scale  C-9  canopy  has  been  used  in  a  variety  of  experiments  and  its  opening 
behavior  was  predicted  with  an  axisymmetric  coupled  model  developed  at  Natick.^®  The  present 
modeling  effort  extends  to  a  cluster  of  three  parachutes. 

The  terminal  descent  shape  from  the  axisymmetric  model  was  used  to  define  the  siuface 
of  the  canopies.  Because  of  the  symmetry  of  this  problem,  a  computational  mesh  was  generated 
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around  one  of  the  three  canopies  to  reduce  the  required  computer  time  and  memoiy  for  the 
numerical  computation s-  An  expanded  view  of  tire  3-D  grid  around  the  parachute  cluster  is 
shown  in  Figure  2,  Also  included  arc  the  shaded  parachute  body  surfaces.  The  full  grid  consists 
of  12  zones  or  blocks.  It  is  extremely  difficult  to  generate  the  grid  for  this  case  using  one  or  a 
very  few'  zones.  I'hc  full  grid  consists  of  a  total  of  84,690  grid  points.  The  grid  points  are 
clustered  near  the  body  surface  for  viscous  turbulent  flow  computations.  The  outer,  in-flow,  and 
down-stream  boundaries  are  placed  sufficiently  aw'ay  from  the  body  surface  that  they  do  not 
interfere  w'ith  the  convergence  and  accuracy  of  the  computed  flow  field  results.  Since  it  is  rather 
time  consuming  to  obtain  a  structured  grid  for  this  problem,  one  has  to  assess  the  use  of  other 
alternate  grid  generation  methods,  such  as  the  unstructured  grid  scheme  and  the  Chimera  overset 
grid  teclmique  for  structured  grids,  Tire  unstructured  grid  generation  taJtes  considerably  less  time 
and  was  assessed  first,  The  surface  grids,  both  structured  and  unstructured,  of  the  parachutes  in 
a  cluster  are  shown  in  Figure  3.  After  the  surface  grids  are  obtained,  the  outer  boundaries  are 
defined,  and  the  3-D  volume  grids  are  generated.  The  grid  generation  process  is  easier  wMt  the 
unstructured  grids.  Therefore,  very  little  effort  would  be  required  to  modify  the  parachute 
cluster  geometry  and  use  the  teclmique  to  predict  cluster  configurations  of  2,3,4  or  more  canopies 
of  any  diameter. 


Figure  2.  Computational  Grid  for  a  Parachute  Cluster. 
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Figures.  Surface  Grids  f  structured  and  unstructured^ 


4.  RESULTS 

The  Input  velocity  used  is  1 8  feet  per  second,  which  was  the  terminal  descent  velocity 
determined  from  the  experimental  results  and  is  very'  elose  to  the  terminal  descent  velocity  used  in 
the  ax i symmetric  models.  All  numerical  computations  wfere  performed  at  this  in-flow  velocity  and 
at  a  =  0“.  Computed  results  have  been  obtained  for  the  3-D  cluster  configuration  and  are  now 
presented.  The  computed  results  were  obtained  using  the  flow  solvers  described  earlier.  Although 
computations  were  performed  on  one  parachute  because  of  symmetry,  the  cluster  consists  of  three 
canopies  wirich  are  shown  in  Figure  3.  Figure  4  show's  the  particle  traces  in  the  vicinit},'  of  the 
parachute.  The  incoming  flow  stagnates  at  the  body  surface,  and  a  separated  flow'  region  is  formed 
in  front.  It  also  shows  How  separating  at  the  skirt  and  then  forming  a  large  region  of  recirculatory 
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flow  behind  the  body.  The  flow  in  the  wake  is  tlu'ee- dimensional  and  shows  the  strong  interactions 
in  the  wake  region,  This  separated  flow  region  in  the  wake  gives  rise  to  lower  surface  pressure  on 
the  outer  body  surface, 


Figure  4.  Particle  Trace  for  Parachute  Cluster,  cc  -  0^. 

The  surface  pressure  distributions  on  the  inner  and  outer  surfaces  of  the  canopies  are 
shovvn  in  Figures  5  and  6,  respectively.  This  case  corresponds  to  a  parachute  polar  angle  of  25  , 
i.e.,  the  canopies  are  25°  away  from  the  vertical  axis.  As  expected,  the  inner  surface  pressure  (see 
Figure  5)  is  quite  uniform  and  is  a  lot  higher  than  the  predicted  pressure  on  the  outer  surface. 

The  pressure  on  the  outside  of  the  canopies  is  low';  also,  there  is  a  variation  in  the  pressure  from 
the  apex  to  the  skirt  of  the  canopies  (see  Figure  6).  The  difference  in  pressure  between  the  outer 
and  the  inner  surfaces  gives  rise  to  drag  force,  which  is  consistent  w'ith  the  payload  for  this 
parachute  at  this  terminal  velocity.  Currently,  the  canopy,  suspension  lines,  and  riser  weights 
are  incorporated  into  the  payload  weight. 
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Figure  5.  Computed  Inner  Surface  Pressure,  a  =  0^.  ( 25^). 


Figure  6.  Computed  Outer  Surl'ace  Pressitre,  a  =  0°.  f  ^  =  25°), 

Figitrc  7  shows  the  pressure  contours  along  tw'o  planes  cutting  through  the  canopies.  For 
the  incoming  flow,  the  pressure  is  uniform  (free  stream  pressure),  As  the  flow  approaches  the 
inner  surface,  it  stagnates  and  forms  a  large  region  of  high  pressure  in  the  vicinity  of  the  inner 


to 


surface.  The  pressure  in  the  wake  region  is  a  low-pressure  region.  Although  there  is  some 
variation  of  pressure  in  the  low-pressure  region,  the  outer  surface  pressure  is  a  lot  lower  than  the 
inner  surface  pressure.  Sharp  changes  in  the  pressure  field  can  be  observed  in  going  from  the 
inner  to  the  outer  surface  because  of  flow  expansion. 


Figure  7,  Computed  Pressure  Contours,  a  =  0°.(^=  25'^y 


The  computed  pressure  map  for  the  parachute  cluster  is  shown  in  Figure  8  for  the  polar 
angle  of  40°,  Again,  it  shows  high  pressure  on  the' inside  surface,  which  is  rather  uniform  except 
for  the  skirt  region.  The  computed  pressure  on  the  outer  surface  is  low^er  and  is  not  as  nniform. 
As  expected,  the  flow  field  is  different  compared  to  the  25°  case.  The  pressui-e  differential 
between  the  inner  and  the  outer  surface  again  determines  the  drag  force  for  the  parachute  cluster 
at  the  prescribed  terminal  descent  velocity, 

As  described  in  an  earlier  section,  the  forces  and  moments  for  the  parachutes  in  cluster  arc 
determined  from  the  computed  solutions.  Figure  9  shows  the  cluster  efficiency  factor  as  a 

function  of  the  polar  angle,  The  cluster  efficiency  factor  is  defined  as  the  ratio  of  the  drag  of 
parachutes  in  a  cluster  to  the  drag  of  a  single  parachute  times  the  number  of  parachutes,  i,e,,  Cpj 

(cluster)  /  [Cpj  (single)  *  number  of  parachutes] .  The  cluster  efficiency  factor  is  usually  less  than 

1 .0  because  of  wake  flow  field  interactions.  Figure  9  clearly  shows  the  computed  cluster 
efficiency  factor  in  the  range  of  0,82  to  0,92,  which  is  similar  to  the  expeomental  data  in  Knacke^ 


for  larger  solid  round  flat  ciiciilai  parachutes.  The  cluster  efficiency  drops  a  little  with  tire 
increase  in  the  polar  angle.  Computed  results  for  ^  ^  25'’  show  that  both  the  structured  and 
unstructured  solutions  produce  almost  the  same  cluster  efficiency  factor. 


Figure  8.  Computed  Pressure  Contours.  Ot  =  0'^.  f  ^  =  40^1. 


Figure  5,  Cluster  Efficiency  Factor  Versus  Polar  Angle'  ^ 
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Figures  10  and  1 1  show  the  computed  force  and  moment  coefficients  as  a  function  of  the 
polar  angle  of  the  clustered  parachute  configuration.  These  results  are  obtained  from  the 
unstructured  flow  computations  for  all  three  parachutes  in  a  cluster.  The  three  force  components 
acting  on  the  parachutes  are  shown  in  Figure  10.  Here,  the  vertical  force  Cy  is  the  drag  force  and 
as  expected,  is  the  largest  component.  The  predicted  drag  is  almost  constant  with  the  polar  angle. 
It  is  seen  to  decrease  only  a  little  with  increased  polar  angle.  The  force  components  in  the  other 
two  directions  are  rather  small.  It  is  interesting  to  note  that  both  force  components  in  the  x  and  z 
directions  are  zero  between  polar  angles  of  30°  and  35°.  Figure  1 1  shows  the  variation  of  the 
moment  components  in  three  directions  with  the  polar  angle.  The  moment  coefficients  are 
referenced  to  the  origin  (0,0,0).  All  three  moment  components  go  through  zero  between  polar 
angles  of  30°  and  35°.  In  particular,  the  component  in  y-direction  is  zero  at  (j>  =  33°.  The 
expected  stable  configuration  for  the  cluster  geometry  is  determined  by  the  condition  at  which  the 
net  moment  about  the  payload  (origin)  is  zero.  In  this  case,  it  happens  at  approximately  at  (j)  = 
32°.  This  computed  polar  angle  compares  very  well  with  =  33°,  which  is  the  experimentally^ 
observed  polar  angle  for  the  stable  configuration.  In  the  experiments,  the  positions  of  the  three 
canopies  were  measured  during  the  opening  process,  and  measured  polar  angle  increased  during 
the  opening,  until  it  stopped  at  around  (j)  =  33°  when  steady  state  was  achieved. 


Figure  10.  Force  Components  Versus  Polar  Angle.  ^ . 
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Figure  11.  Moment  Components  Versus  Polar  AnaLe..^  ^ . 


5.  CONCLUDING  REMARKS 

A  computational  study  has  been  performed  to  determine  the  aerodynamics  of  a  cluster  of 
round  parachutes  using  CFD.  The  results  given  here  are  the  first  predictions  of  descent 
characteristics  for  a  cluster  of  three  half-scale  C-9  parachutes.  In  particular,  the  results  include 
the  aerodynamic  flow  field  and  geometry  of  the  parachute  cluster  with  assumed  shapes  for  the 
individual  canopies.  Using  periodic  symmetry  of  the  cluster  configuration,  the  pressure  over  the 
inner  and  outer  surfaces  Lf  a  single  canopy  in  the  cluster  is  used  to  calculate  the  net  vertical  and 
radial  forces  acting  on  the  canopy.  A  manual  iterative  process  is  used  to  determine  the  expected 
stable  configuration  for  the  cluster  geometry  by  determining  the  condition  at  which  the  net 
moment  about  the  payload  (origin)  is  zero.  The  corresponding  angular  location  of  the  canopies 
predicted  by  the  CFD  cdmputations  is  compared  with  the  available  experimental  results  and  is 
found  to  agree  well  with  the  data.  It  is  shown  that  significant  progress  has  been  made  in 
determining  the  terminal  descent  flow  field  characteristics  of  a  particular  parachute  cluster 
configuration.  These  ccjmputational  solutions  provide,  for  the  first  time,  an  understanding  of  the 
flow  field  in  and  around  parachute  clusters. 

This  report  has  presented  the  current  status  of  the  modeling  effort  to  address  the 
complexity  of  the  parachute  characteristics  and  the  associated  3-D  flow  fields  of  clusters  of 
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round  parachutes.  Future  research  will  include  the  use  of  Chimera  overset  grid  approach  to 
compute  the  time-dependent  fluid-dynamic  forces  associated  with  “rigid”  parachutes  in  a  cluster. 
This  technique  will  allow  the  determination  of  the  steady  state  stable  configuration  for  a  given 
cluster  design  “dynamically”  within  one  numerical  simulation.  The  developed  technique  is 
expected  to  assist  in  the  development  of  future  U.S.  Army  airdrop  systems  and  other  roimd 
parachute  systems.  The  capability  of  accurately  predicting  the  behavior  of  parachute  systems 
will  significantly  reduce  the  amount  of  testing  currently  required. 
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collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  121 5  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 
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